A study of the anisotropy associated with dipole asymmetry in heavy ion colhsions 
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The anisotropy associated with the initial dipole asymmetry in heavy ion collisions is studied via 
the two-particle relative azimuthal azimuthal angle {Acj> — (jf — (t>°) correlations, within a multi-phase 
transport model. For a broad selection of centrality, transverse momenta {p'^ ) and pseudorapid- 
ity {rf^'^), a fitting method is used to decompose the anisotropy into a rapidity-even component, 
characterized by the Fourier coefficient Vi, and a global momentum conservation component. The 
extracted ni values are negative for pt ^ 0.7 — 0.9 GeV, reach a maximum at 2-3 GeV, and de- 
creases at higher pt- The vi values vary weakly with tj and centrality, but increases with collision 
energy and parton cross-section. The extracted global momentum conservation component is found 
to depend on Ar; = rf" — rf' for |Ar;| < 3. 

PACS numbers: 25.75.-q 



10 A large azimuthal anisotropy is observed for parti- 

11 cle production in heavy ion collisions at the Relativis- 

12 tic Heavy Ion Collider (RHIC) and the Large Hadron 

13 Collider (LHC). This anisotropy is understood as the re- 

14 suit of pressure-driven collective flow, which converts the 

15 spatial asymmetries of the initially created matter into 

16 transverse momentum (px) space [l|. Such anisotropy 

17 can be expressed as Fourier series in azimuthal angle (0): 



dNS 
~d4 



1 + 2 ^ w„(pt) cosn((/i - $„), 



(1) 



n=l 



18 where Vn and $„ represent the magnitude and direction 

19 of n*'^-order flow respectively. Asymmetries giving rise to 

20 non-zero t;„ are associated with either an average shape 

21 (for n = 2) or shapes due to the spatial fluctuations of 

22 the participating nucleons [2t5| ; the n = 1 term is known 

23 as dipole asymmetry. Due to the approximate boost in- 

24 variance of the created matter, the u„ are expected to 

25 be even functions in rapidity (rapidity-even) [6|, |7[ . How- 

26 ever, a small rapidity-odd u„ component, especially for 

27 n — 1, may arise from the "sideward" deflection of the 

28 created matter [S\ . This rapidity-odd vi is expected to be 

29 very small for the rapidity range and the collision ener- 



111, and is not considered in this 



30 gies considered here 

31 paper. 

32 The harmonics in Eq. [T] is also reflected in the two- 

33 particle correlation (2PC) in A(f) ^ (p^ — (f>^: 



dN, 



pairs 



dAcj) 



l + 2j2vn,niPT,PT)cOSnA(j) , (2) 



n=l 



with 



Vn.niPT^PT) = W„(pt)w„(Pt) 



(3) 



Recent studies at the LHC show that the latter factor- 
ization is valid for n > 2 at low px for particle pairs 
with a large rapidity gap 12, ll3|, but fails for n = 1 
over the entire pT range. This is understood to be due 



39 to global momentum conservation for a finite multiplic- 

40 ity system, which induces a negative dipole component 

41 inEq.|3][i,[iM3|: 



vi,i{pt,Pt) = vi{Pt)Mpt) 
1 



cPtPt 



KM{p\) 



(4) 
(5) 



42 where M is the total multiplicity in the event and K is 

43 the fraction of particles that are correlated via momen- 

44 tum conservation. A recent study by the ATLAS CoUab- 

45 oration |l3[ shows that wi^i can indeed be described by 

46 Eq. m over broad ranges of p^ . The vi values extracted 

47 via a two-component fit are negative for px ^ 1 GeV, in- 

48 crease to a maximum at 4-5 GeV, and decrease at higher 

49 px- A similar fitting method was previously applied by 

50 Retinskaya et.al ^ to the ALICE 2PC data [1^], al- 

51 though the data have less precision and may not be free 

52 of short range correlations due to smaller A77 gap. 

53 Many models based on viscous hydrodynamics or 

54 transport theory have been developed to describe the u„ 

55 and Vn,n data [6, 20 27]. They show that sizable higher- 

56 oder flow {n > 2) can be generated by including fluctu- 

57 ations of the initial geometry. In some cases, reasonable 

58 agreement with the Vn data has been achieved. However, 

59 due to the complication associated with the global mo- 

60 mentum conservation, most of these model attempts have 

61 avoided a prediction for vi. Recently, Gardim et.al [6| 

62 performed the first calculation of the rapidity-even f 1 in 

63 a hydrodynamic code by directly cancelling the momen- 

64 tum conservation from ui^i by chosing appropriate weight 

65 for each particle (the cancellation is exact if the global 

66 conservation is linear in px); a modified event plane (EP) 

67 method based on this cancellation procedure was also 

68 proposed [17| . Note that the systems that are related via 

69 momentum conservation could be smaller than the whole 

70 event (for example a dijet system) [28|. Therefore, the 

71 value of c (Eq. [SJ may depend on other kinematic vari- 

72 ables such as the A77 of the particle pairs. Hence, it is 
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FIG. 1: (Color online) 
1'i,i(Pt.Pt) (symbols) 
from AMPT (left pan- 
els) and HIJING (right 
panels) for b = 4 f m (top 
panels) and 6 = 8 fm 
in AuAu collisions at 
.JsWn = 0.2 TeV. Lines 
in the left (right) panels 
are result of a global fit 
to Eq. |4] (a fit to pure 
momentum conservation 
component: cpxPt)- 



advantageous to determine c via the data driven method 
used by the ATLAS Collaboration. 

In this work, we present a quantitative study of the 
rapidity-even v\ using a generalization of the ATLAS fit- 
ting method within the AMPT (A Multi Phase Trans- 
port) model [29|. This study allows us to make quantita- 
tive predictions on how v\ varies with px, '^^ centrality, 
collision energy and partonic cross-section in a transport 
framework, as well as to compare with the ATLAS data. 

The AMPT model has been successfully used to under- 
stand the higher-order Vn associated with initial asym- 



metries [20|, l2 l |. | 2 3l|. It combines the initial condition 



from HIJING [30| and final state interaction via a par- 
ton and hadron transport model. The initial condi- 
tion from HIJING is seeded by a Glauber model for 
nucleon-nucleon scatterings, hence it contains the fluctu- 
ations which generate the event- by- event spatial asym- 
metries. The nucleon-nucleon scatterings provide both 
minijet partons for hard scattering and strings for soft co- 
herent interactions; these strings are then converted into 
soft partons via a string melting scheme. The parton and 
hadron transport is responsible for transforming the ini- 
tial asymmetries into the momentum anisotropy. The 
parton transport include only elastic scattering whose 
cross-sections are controlled by the values of the strong 
coupling constant (us") and the Debye screening mass (/i): 

^™' T « %WJ{2^?) (6) 



do 



2(t-Ai2)2 



where t is the standard Mandelstam variable for four mo- 
mentum transfer. Note that the total cross-section can be 
increased either with a larger as or a smaller value of fi. 
However changing as is more effective for momentum dis- 
sipation than changing /i, as decreasing fi includes only 



softer scattering (Eq. [5]). Partons are recombined into 
hadrons at the freezeout, followed by a hadronic trans- 
port. 

AMPT events are generated for two different as values 
for parton transport: set-A with as — 0.47 and set-B 
with as — 0.33 from 2lj 3l|. We vary the total cross- 



section from 1.5 mb to 10 mb by adjusting /z, separately 
for set-A and set-B. Varying as and /z allows us to study 
how vi depends on the strength of final state interactions. 

13 These studies are carried out at RHIC energy (Au-|-Au 

14 at ^s]^ = 0.2 TeV) and/or LHC energy (Pb+Pb at 

15 \/snn = 2.76 TeV) for the impact parameters 6 = 4 fm 

16 and b = 8 fm, corresponding to approximately 0-10% and 

17 30-40% most central events, respectively. 

18 Wc follow the procedure of extracting vi as docu- 
mented in 13]. First, the correlation function for a given 

20 p^, pJ^ and A?/ range is constructed from all final state 

21 particles, and vi,i{p^, p!^, Arj) values are calculated. Sec- 

22 ond, a least-square minimization of these wi^i values fol- 

23 lowing Eq. 2] is used to obtain uf'*(pT) for a given Ary 

24 range, where the wf'*(pT) function is defined by its val- 

25 ues at 9 px points (0.5,0.7,0.9,1.1,1.5,2.0,3.0,3.8,4.6 GeV) 

26 via a cubic-spline interpolation procedure. 

27 Figure [1] (a) and (c) show examples of such a study 

28 at RHIC energy, with the fit overlaid with the wi.i from 

29 AMPT for set-A with a lOmb partonic cross-section. A 

30 rapidity gap of 1.5 is imposed to suppress the short range 

31 correlations. The ui i values for different p^ selections 
cross each other and change sign as one increase p^ , a 
feature qualitatively similar to the ATLAS data. 

Figure[T](b) and (d) show the 111,1 values obtained from 
HIJING, which has no final state interaction. They vary 
nearly linearly with p^ x pJj., consistent with behavior 
expected for global momentum conservation. This sug- 
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FIG. 2: (Color online) uf " vs pt at RHIC energy and LHC. 



138 gests that the complex dependence of the v\;\ on 'p^ 

139 in the left panels is a natural consequence of final state 

140 interactions. 

141 Figure [2 shows wf'*(pT) extracted from the fit. The 

142 t)f'' function crosses zero at low px, with the crossing 

143 point of pt,o = 0.7 — 0.8 GeV at RHIC energy and 0.9 

144 GeV at LHC energy, respectively. Since the integral 

145 V\ should be zero: (pt'^iIpt)) — 0, the change of pt,o 

146 mainly reflects the change of the (pt) and (p^) (value 

147 shown in the legend). One can show that if z;i(pt) is a lin- 

148 ear function of px, the crossing point is pT,o — (Pt)/(pt) ■ 

149 However, a change in the (px) and (p^) should have lit- 

150 tie influence on the v\ values at high px; there the v\ 

151 values are more sensitive to transport properties of the 

152 medium. The wf''(pT) increases with px to a value ^0.1 

153 at 2-3 GeV, and then decreases at higher px. The wf'* 

154 values for 6 = 4 f m and 6 = 8 f m are comparable, consis- 

155 tent with the lack of centrality dependence in the initial 

156 dipole asymmetry [J, |5| . The ?;f '* value at LHC energy is 

157 larger, consistent with a stronger collective flow at LHC. 
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FIG. 3: (Color online) wf'* vs px for set-A (left panels) and 
set-B (right panels) for impact parameter 6 = 8 fm, and com- 
pared with ATLAS Pb-|-Pb data. 



rj . The resulting 4-D dataset vi,i{p^,p^,r]^,ri ) is then 
parameterized by: 

/a b a b\ Fit/ a a\ Fit/ b b\ 

Vi,i[Pt,Pt,V ,?? ) = Wi {Pt,V m {Pt,V ) 



f a b\ a b 

'C(j] ,V )PtPt 



(7) 



159 Figure [3] shows ff''(px) for set-A and set-B at different 

160 parton cross-sections tuned by changing /z. The mag- 

161 nitude of wf''(px) clearly grows with increasing cross- 

162 section. However the vi values are more sensitive to 

163 changes in as'- a change of a^ — 0.47 for set-A to 0.33 for 

164 set-B almost leads to a 60% increase of the peak value of 

165 vi. Results for set-B with a 1.5 mb cross-section agree 

166 better with the ATLAS data at low px. However both 

167 sets fail to describe the trend of the data for px > 3 

168 GeV. This reflects a limitation in the AMPT model for 

169 predicting the behavior of u„ at high px: the jet quench- 

170 ing physics is important at high px, but is not modelled 

171 correctly in AMPT. Similar studies for V2 also found that 

172 the set-B 1.5 mb has the best agreement with the data 

173 forpx <2 GeV [21|. 

174 The original fitting method of ATLAS can be gener- 

175 alized to measure the rj dependence of vi in AMPT. We 

176 first perform the 2PC analysis in bins of p^,p!j., 77*^, and 



to directly obtain wf'*(px,?y). A minimum jAryj separa- 
tion for the pairs is imposed to suppress the short range 
non-fiow correlations. Therefore the total number of wi.i 
data points used in the fit depends on this gap require- 
ment and the binning in px and 77. 

In the default fit, both t^f'* and c functions are as- 



sumed to be symmetric in 77 i.e. 



,Fit 



iv) 



,,Fit 



i-v) 



and c{rj^,ri^) = c{—rj'^,—T]^). In addition, we also have 
c{r]'^,r]°) = c{rj^,7f). We divide the |7;| < 3 into 12 
equally spaced rj bins and use the same 9 bins in px 
as for Fig [1] and O The fit uses al2x9=108 mesh 
in (7/,px) space with a simple linear interpolation, which 

91 has 54 independent mesh nodes due to symmetry in j]. 

92 This together with 24 parameters for c, gives a total of 

93 78 fitting parameters. The quality of the fit is generally 
comparable to those shown in Fig. [1] 

Figure [4] shows the rj dependence of the vi obtained 
this way for two px bins. The error bands account for 
the statistical uncertainty as well as the variation of the 
fit result by changing the |A?7| gap from 0.5 to 2.5. The 

199 vf '*(??) generally shows a weak ry-dependence at RHIC 

200 energy, especially in the more central collisions. At LHC 

201 energy, however, the wf '*(??) shows a small dip at mid- 
202 rapidity. The physics behind this dip is currently under 

203 investigation, but probably is related to the stronger lon- 

204 gitudinal flow at higher collision energy. Figure [5] shows 

205 the values of c^q^.if') plotted as a function of IA77I for 

206 6 = 8 fm at RHIC and LHC. They fall approximately on 

207 a common curve for |A?7| > 2, suggesting that c nearly 

208 depends only on |A?7| for pairs with large rapidity sep- 

209 aration. The values of c at large JAt^J are not sensitive 

210 to the IA77I gap used in the global flt (compare the sohd 

211 and open symbols), and they are three times larger than 

212 that calculated via Eq. [5] assuming K = I at LHC (right 

213 panel), but are comparable at RHIC (left panel). 
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15-2 GeV -"^ by the hydrodynamic models. The vi function shows 
-)43 very httle centrality dependence but increases with colli- 
n M4 sion energy and parton cross-sections, suggesting that its 
us magnitude increases with the strength of the final state 
interactions. By choosing the parameters in AMPT care- 
er fully, reasonable agreement with the ATLAS data can be 
MS achieved for px < 2 GeV, but not at higher px- The 
.M9 fitting method is extended to directly extract the rj de- 
-250 pendence of vi and the global momentum conservation 

251 component (Eq. [7]). The extracted vi shows a weak de- 

252 pendence on rj. The coefficient of global momentum con- 

253 servation component is found to depend on Ary and differs 

254 from that given by Eq. [5l This extended method should 

255 allow the extraction of the vi (pt , v) ^^'^ validation of the 

256 EP method (It requires detector acceptance to px = 0, 

257 which is usually not satisfied in actual experiment). 

258 We acknowledge valuable discussions with Matthew 

259 Luzum and Jean- Yves OUitrault, in particular, the re- 

260 lation between the event plane method and two-particle 

261 correlation method. We thank Roy Lacey for a careful 

262 proofreading of the manuscript. This research is sup- 

263 ported by NSF under award number PHY-1019387. 



FIG. 5: c{r)'',r]°) vs. jry^ - r)^\ for A77 > 0.5 and A77 > 2 
obtained from the Global fit. Data points at given |A?7| value 
corresponds to different rf and r]" combinations 



214 In the original functional Eq.|4l parameter c is assumed 

215 to be a single number for each centrality class. Figure [5] 

216 clearly show that this is not the case in AMPT, and 

217 the global momentum conservation contribution changes 

218 with IA77I. Some of this dependence is due to the influ- 

219 ence of the short range correlations such as jet fragmen- 

220 tation and resonance decay, which tend to give a negative 

221 contribution to c at small IA77I. However, the values of c 

222 still change with |A?7| up to jATyj = 3 then seem to flat- 

223 ten out at higher |A7y|, but at a level much bigger that 

224 predicted by Eq. [5l This may suggest that evaluating c 

225 with all particles in the event might be incorrect. In the 

226 modified EP method of Ref. [17|] , a px dependent weight 

227 ui = px — (Px)/(pt) is calculated for each particle used 

228 in the EP calculation. Note that if c depends on 77 of the 

229 two particles, {p\) and (px) need to be calculated in bins 

230 of 77, instead of a single number calculated over the full 77 

231 range of the EP. In other words, the weight is a function 

232 of px and 77: w(j)T,r]) = px - (Pt) iv) / (pt) iv) ■ 

233 In summary, we study the vi associated with the ini- 

234 tial dipole asymmetry in the AMPT Monte Carlo model 

235 using wi,i = (cos A0) obtained from the two-particle cor- 

236 relation technique. A global fit of the wi^i using the pro- 

237 cedure from the ATLAS Collaboration 'IS"] allows the 

238 simultaneous extraction of 7Ji(px) and the global momen- 

239 tum conservation component. The extracted vi function 

240 is negative at low px, and increases with px until 2-3 

241 GeV, qualitatively agreeing with the behavior predicted 
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